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Using a dynamic cluster quantum Monte Carlo approximation, we study a two-dimensional Hub- 
bard model with a small orthorhombic distortion in the nearest neighbor hopping integrals. We 
find a large nematic response in the low-frequency single-particle scattering rate which develops 
with decreasing temperature and doping as the pseudogap region is entered. At the same time, 
the d-wave superconducting gap function develops an s-wave component and its amplitude becomes 
anisotropic. The strength of the pairing correlations, however, is found to be unaffected by the 
strong anisotropy, indicating that d-wave superconductivity can coexist with strong nematicity in 
the system. 

PACS numbers: 74.72.-h, 74.20.-z, 71.10.-w 



Introduction- Nematic correlations, where the C4 ro- 
tational symmetry of the square lattice is spontaneously 
broken and reduced to C2, have been observed in the 
electronic structure in a number of strongly correlated 
electron materials, including the strontium-ruthenate ox- 
ide compounds, the cuprates and the iron-pnictide ma- 
terials [1]. Despite the presence of very small structural 
anisotropics in the plane, these systems display very large 
nematic anisotropics in properties such as the dc and in- 
frared resistivities 0, the magnetoresistance the 
Nernst signal the neutron scattering spectrum ^ 
as well as in scanning tunneling spectroscopy [gI, . 

From a theoretical perspective, the idea that the pseu- 
dogap in the cuprates may be associated with a ne- 
matic phase, which arises from melting of a stripe or- 
dered phase, was first discussed in Ref. [8] and later 
expanded on in Ref. [9]. A first microscopic theory 
of how a nematic phase can emerge from an isotropic 
Fermi liquid thr oug h a Fermi surface instability was pre- 
sented in Ref. [10|. And indeed, both strong-coupling 
RVB slave-boson mean-field calculations for the 2D t-J 
model [11, 12] and weak-coupling functional renormaliza- 
tion group studies jlsl-flBj of the 2D Hubbard model have 
found a Pomeranchuk instability towards a d-wave Fermi 
surface deformation (for a review see Ref. [16]). This ne- 
matic instability has been found to compete with d-wave 
superconductivity, a picture that has also emerged from 
exact diagonalization [l7| and variational Monte Carlo 
studies [18| of the t-J model. 

Here we are interested in the question of how a small 
orthorhombic distortion that exists in the system affects 



the electronic structure in a strongly correlated electron 
system. Consistent with experiments, recent theoretical 
studies have found a large nematic response to a small 
structural anisotropy. In particular, a recent study of the 
2D Hubbard model using cellular dynamical mean-field 
theory and a dynamic cluster approximation (DCA) of a 
small 2 X 2-site cluster has found a large nematic response 
in the low-frequency behavior of the single-particle self- 
energy as well as in the conductivity when the system is 
doped towards the Mott insulating state [22]. This study 
was restricted to zero temperature and a small 2x2 clus- 
ter. Here, using a larger 4x4 cluster DCA calculation of 
a similar 2D Hubbard model with a small anisotropy in 
the single-particle hopping integral, we study the tem- 
perature dependence of the nematic response as well as 
its connection with the pseudogap and its interplay with 
the superconducting behavior of this model. 



We find a large nematic anisotropy in the single- 
particle scattering rate, which increases with decreasing 
temperature and doping. We show that this nematic 
anisotropy becomes significant inside the pseudogap re- 
gion below a temperature T*(^), pointing to a deep re- 
lationship between the pseudogap and nematic correla- 
tions. Moreover, we find that the orthorhombicity in the 
single-particle hopping leads to an anisotropic supercon- 
ducting gap, with both d-wave and s-wave components 
similar to what was discussed in Ref. Q , but has no effect 
on the superconducting transition temperature. 
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Formalism- The orthorhombic 2D Hubbard model we 
analyze [22], 

has a nearest neighbor transfer integral tij = —t^/y for 
hopping along the x- and ^/-direction, respectively, a 
Coulomb repulsion U and chemical potential /i. As in 
Ref. [ill, we use a small anisotropy t^/y = t(l ± ^/2) 
with ^ = 0.04 to break the C4 lattice symmetry of the 
model. In the following we will measure energies in units 
of t and set the Coulomb repulsion U = 6t unless other- 
wise noted. 

We will use the DC A to study the model in Eq. ([1]). 
The DCA [l9|, [20| maps the bulk lattice problem onto 
an effective periodic cluster embedded in a self-consistent 
dynamic mean-field that is designed to represent the re- 
maining degrees of freedom. The DCA calculations were 
carried out on a 4x4 cluster and the effective cluster 
problem was solved using a Hirsch-Fye quantum Monte 
Carlo algorithm [21 1. 

Results- In contrast to the 2 x 2-site cluster study per- 
formed in [22], we have found that, in the 4x4-site 
cluster, an additional next-nearest neighbor hopping t' 
reduces the nematic response relative to that found for 
f = 0. We therefore set = in the 4x4-site cluster 
calculations. 

In order to examine how the directional anisotropy in 
the hopping integral affects the Fermi surface of the cor- 
related system, we plot in Fig. 1 the gradient of the 
momentum distribution function, |Vn(k)| with n(k) = 
^(^iccr^kcr)' ^^^^ quadrant of the Brillouin zone. We 

have interpolated the DCA cluster self-energy S(K, iun), 
which is defined on a 4 x 4 grid of K-points using a smooth 
spline. 

The panels on the left show the results for the 4x4- 
site cluster for U = 6t and t' = 0. Panels a) through 
c) show results for different dopings 5 calculated at a 
fixed temperature T = 0.074t, while panel d) shows the 
result for a lower temperature T = 0.05t and a doping of 
S = 0.05. From this one sees that at high temperatures 
and doping levels, |Vn(k)| is fairly symmetric. But at low 
temperatures and dopings, directional anisotropy clearly 
develops in both its shape and, much more noticable, 
in its magnitude, with a larger value in the re gion near 



4x4 cluster 
U=6t, t' = 



(5=0.05 



2x2 cluster, T=0.074t 



(0,7r) than near (7r,0). As pointed out in Ref. [22|, this 
refiects a stronger anisotropy in the imaginary part of the 
self-energy than in the real part. 

This result is supported by the plots of the real and 
imaginary parts of the self-energy in Fig. 2. There one 
clearly sees that the imaginary part exhibits qualitatively 
different behavior between (tt, 0) and (0,7r), while the 
real part shows much less anisotropic. For K = (0,7r), 
|ImS(K,a;n)| displays Fermi liquid like behavior with a 
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FIG. 1: (Color online) Plots of the gradient of the k-space 
occupation, |Vnfc| in the first quadrant of the Brillouin zone, 
calculated for the system with small directional anisotropy. 
On the left, results are shown for the 4x4 cluster with U — 6t 
and = in panels a) through c) for a fixed temperature 
T = 0.074t but varying doping as indicated, and in panel 
d) for a lower temperature T = 0.05t and a doping S = 0.05. 
The right side illustrates the effects of cluster size and varying 
parameters for a doping S = 0.05. Panels e) and f) show 
results for a 2x2 cluster for different values of U and t' as 
indicated. Panels g) and f) show results for the 4x4 cluster 
for t' = for different values of U. 



downturn at low frequencies, while for K = (7r,0), it 
keeps rising with decreasing cun- A similarly strong angle 
dependence of the self-energy in the nematic phase was 
also found in Ref. [10]. In addition, one sees that the 
anisotropy in the imaginary part goes away at high fre- 
quencies, consistent with the behavior found in Ref. [22|. 

We generally find that the nematic response in the 
4x4-site cluster is not larger than that found in the 2x2- 
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FIG. 2: (Color online) Real and imaginary parts of the self- 
energy I1(K, ujn) versus frequency uJn for K = (tt, 0) and (0, tt) 
for a doping 5 — 0.05 and temperature T — 0.125t. Results 
for both the anisotropic and isotropic cases are shown. 



site cluster. As an example, we show on the right hand 
side of Fig. 1 in panels e) and f) results for the 2x2-site 
cluster for U = and = and /7 = lOt and t' = -0.3t, 
respectively. Here one sees that in the 4-site cluster, the 
nematic response is dramatically increased for larger val- 
ues of U and a finite = —0.3t. In the 16-site cluster, 
however, this does not seem to be the case. While a sys- 
tematic study of larger values of U and finite f is not 
possible in the 16-site cluster, because of the QMC sign 
problem, we show as an example in panels g) and h) 16- 
site cluster results for two different values of /7, U = 6t 
and 8t, respectively. There is no obvious enhancement of 
the response for the larger value of U. 

In order to further examine the temperature and dop- 
ing dependence of the nematic response, we have calcu- 
lated the quasiparticle weight approximated by 



Im 



S(K,7rT) 



ttT 



(2) 



In Fig. 3 we plot the nematic anisotropy of Zk defined 
as 



crz 



^(0,7r) 



^(7r,0) 



^(0,7r) 



^(7r,0) 



(3) 



for the 4 X 4-site cluster. One clearly sees that the nematic 
response increases with decreasing temperature and dop- 
ing, reaching remarkable levels of over 50% at the lowest 
temperature and doping level we have studied. 

Also indicated in this plot by the solid blue line is the 
pseudogap temperature T*(^) as in [23]. Below this tem- 
perature, a pseudogap opens in the single-particle den- 
sity of states as well as in the spin excitations. Here we 
have determined T* {S) from the temperature T at which 
the uniform magnetic susceptibility Xs(q = 0?^) has a 



maximum 



24| . It is interesting to note that this temper- 
ature is not affected by the orthorhombic anisotropy in 
the Hamiltonian. Conversely, as evidenced by this plot, 
the emergence of a strong nematic response is intimately 
linked to the opening of the pseudogap. Outside the 
pseudogap region, the nematic anisotropy is small. With 
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FIG. 3: (Color online) Color map of the anisotropic ratio 
of the quasiparticle weight, az, over the temperature-doping 
plane, for U=6t. The blue curve indicates the pseudogap 
temperature T*(6) which is obtained as the temperature at 
which the uniform magnetic susceptibility Xm{q = (0,0), T) 
has a maximum. 



decreasing temperature and doping, it increases dramat- 
ically along the pseudogap temperature T*(^) and then 
saturates deep inside the pseudogap region. This points 
towards the pseudogap being a necessary precondition 
for nematic physics to appear in this system. 

Finally, we note that we do not find spontaneous sym- 
metry breaking to a nematic state in the C4 symmetric 
model, i.e. for ^ = 0, for the temperatures we have stud- 
ied. Note, however, that this does not exclude the pos- 
sibility of an instability at a lower temperature. In fact, 
the large nematic response we find at the temperatures 
we have studied is strongly suggestive of an instability 
at lower temperatures. A detailed study of this issue is 
beyond the scope of the present paper and will be the 
subject of future work. 

It is interesting to ask how the strong nematic 
anisotropy that is present in the single-particle scattering 
rate affects the superconducting behavior of the system. 
In Ref. [8| it was pointed out that, in a nematic phase, the 
superconducting order parameter will have mixed (1^2 _y2 
and extended s symmetry. This question was also ad- 
dressed in a previous slave-boson mean-field study of an 
anisotropic t-J model in Ref. [ll^. 

Here we study the superconducting behavior of the 
anisotropic Hubbard model ([T]) by solving the Bethe- 
Salpeter equation in the particle-particle channel 



'{K, K')x7{K')MK') = KMK) , (4) 
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FIG. 4: (Color online) (a) Leading d-wave eigenvalues of 
the Bethe-Salpeter equation (|4j) in the particle-particle chan- 
nel versus temperature. The eigenvalues for both the isotropic 
and anisotropic models are plotted, (b) The momentum de- 
pendence of the leading eigenvector (j)d{K^nT) at the low- 
est Matsubara frequency along the path in Brillouin zone 
shown in the inset, (c) The relative s-wave contribution to 
the d-wave contribution, |($d(7r,0) + $d(0, 7r))/($d(7r, 0) - 
$d(0,7r))|, as a function of temperature for different dopings. 



where K = (K,cJn) and Tpp{K,K') is the irreducible 
particle-particle vertex with center of mass momentum 
and frequency Q = 0. The coarse-grained bare particle- 
particle Green's function Xo^{K) = A^c/A^ Ek ^t(K' + 
k.\ujn)Gi{—'K — k, — cjn) is calculated from the lattice 
Green's function Gcr(k, cj^) = [^^n — ^k+M — 5](K, 0;^)]""*^. 
The system undergoes a superconducting transition at 
the temperature Tc where the leading eigenvalue be- 
comes one, and the momentum and frequency depen- 
dence of the corresponding eigenvector (j)a{K) deter- 
mines the symmetry of the superconducting state. At 
Tc, (j)a{K) becomes identical to the superconducting gap 
function. In the isotropic model, the instability occurs in 
the spin singlet, even frequency channel and the leading 
eigenvector has d-wave symmetry [25]. Here we focus on 
the low-doping region, where the nematic response of the 
anisotropic model is most pronounced. 

Fig. 4a shows the temperature dependence of the lead- 
ing eigenvalues of the anisotropic and isotropic models for 
two different dopings 5 = 0.10 and 0.05. In the temper- 
ature range we could access — the QMC Fermion sign 
problem prevents us from going to lower temperatures 
— the leading eigenvalues of the isotropic = 0) and 
anisotropic = 0.04) models are essentially identical. 
This indicates that the strong nematic anisotropy which 
is present in the single-particle scattering rate has essen- 



tially no effect on the strength of the pairing correlations. 
We also calculated the leading eigenvalue in an 8-site 
cluster, which allows us to reach the temperature where 
the eigenvalue crosses one. There we find that the super- 
conducting transition temperature Tc of the anisotropic 
model is unchanged from the isotropic case. 

In Fig. 4b we plot the momentum dependence of the 
leading eigenvector ^(K, ttT) along the line indicated in 
the inset. From this one sees that in the anisotropic 
case, it has predominantly d-wave structure, just as in the 
isotropic model. There is an additional s-wave contribu- 
tion in the anisotropic model, which causes an anisotropy 
between (tt, 0) and (0,7r), but leaves the nodal point un- 
changed. 

In Fig. 4c, we plot the s-wave contribution related 
to d-wave |0((7r, 0), ttT) + ^((0, tt), ttT) |/|0((7r, 0), ttT) - 
0((O, tt), 7rT)| as a function of temperature for different 
dopings. Once again, we find that the anisotropy in the 
eigenvector develops and increases when the system is 
cooled and doped into the pseudogap region. 

These results demonstrate that superconductivity can 
coexist with nematicity in the system. In addition, they 
indicate that the strength of the pairing correlations 
and Tc are essentially unaffected by the large anisotropy 
present in the single-particle scattering rate. This is sim- 
ilar to what was found in earlier slave-boson mean-field 
calculations [llj . 

Conclusions We have used a dynamic cluster quan- 
tum Monte Carlo approximation to study a two- 
dimensional Hubbard model with a small (4%) or- 
thorhombic distortion in the nearest-neighbor hopping 
integrals tx 7^ ty We have found a large nematic re- 
sponse, i.e. a difference in the electronic properties be- 
tween K = (7r,0) and (0,7r), with levels up to 55% in 
the scattering rate Im i;(K, cj^), which develops as one 
enters the pseudogap phase and increases as the doping 
and temperature is further reduced. Similar behavior is 
found in the d-wave superconducting gap function, which 
develops an s-wave component in the pseudogap phase 
which reaches levels up to 40%. In contrast, we find that 
the strength of the pairing correlations is essentially unaf- 
fected by these anisotropics. These results demonstrate 
a close link between the pseudogap and nematic corre- 
lations, and indicate that d-wave superconductivity can 
coexist with a strong nematicity in the system. 
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